Active noise cancellation in a suspended interferometer 

Jennifer C. Driggers/ Matthew Evans, ^ Keenan Pepper,^ and Rana Adhikari-^ 

LIGO Laboratory, California Institute of Technology, Pasadena, CA 9112^^ 
LIGO Laboratory, Massachusetts Institute of Technology, Cambridge, MA 02139 
University of Galifornia, Berkeley, Berkeley, GA 94720 

(Dated: 13 December 2011) 

We demonstrate feed- forward vibration isolation on a suspended Fabry-Perot interferometer using Wiener 
filtering and a variant of the common Least Mean Square (LMS) adaptive filter algorithm. We compare the 
experimental results with theoretical estimates of the cancellation efficiency. Using data from the recent LIGO 
Science Run, we also estimate the impact of this technique on full scale gravitational wave interferometers. 
In the future, we expect to use this technique to also remove acoustic, magnetic, and gravitational noise 
perturbations from the LIGO interferometers. This noise cancellation technique is simple enough to imple- 
ment in standard laboratory environments and can be used to improve SNR for a variety of high precision 
experiments. 
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I. INTRODUCTION 

The next generation of interferometers for 
gravitational-wave detection, including the Laser 
Interferometer Gravitational Wave Observatory (LIGO), 
will have unprecedented sensitivity to astrophysical 
event At low frequencies (~10 Hz) it is likely that 
the displacement noise of the suspended mirrors will 
be limited at the 10~^^ m/>/Hz level by fluctuations 
in the Newtonian gravitational forces^ ^. The sources 
of the fluctuations are density perturbations in the 
environment (e.g. seismic and acoustic) and mechanical 
vibrations of the nearby experimental apparatus. While 
care will be taken to mitigate the sources of all of these 
fluctuations, further reductions of this Newtonian noise 
may be made by carefully measuring the source terms 
and subtracting them from the data stream (offline) 
or in hardware by applying cancellation forces to the 
mirror. 

To demonstrate the efficacy of this technique, we 
demonstrate below that offline subtraction of seismic 
noise can be done using static Wiener filtering based on 
an array of seismic sensors. The online adaptive sub- 
traction is also shown to approach the 'optimal' Wiener 
limit^. Through nonlinear processes, seismic noise be- 
low 10 Hz has been shown to limit the performance of 
gravitational-wave detectors. This technique will prove 
to be of substantial value in reducing the non-stationarity 
of the detectors. 



II. STATIC WIENER FILTERING 

To find a linear filter that will improve a chosen signal, 
we must first define what it means to 'improve' the signal. 
The figure of merit (^) that we use for calculating the 



Wiener filters in this case is the expectation value of the 
square of the error signal (e), where the error signal is 
defined as the difference between the target signal to be 
minimized and the estimate of that signal calculated from 
the filtered witness channels. 
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Here, E[^] indicates the expectation value of *, rep- 
resents the tap weights of the filter, d{n) is the target 
signal which we would like to minimize, p is the cross- 
correlation vector between the witness and target signals, 
and R is the autocorrelation matrix for the witness chan- 
nels. When we solve Eq. [l]by setting 
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Eq.|3]finds the FIR (Finite Impulse Response) filter co- 
efficients which minimize the RMS of the error e by opti- 
mizing the estimate of the transfer function between the 
witness sensors and the target signal. Since the matrix 
R is Block Toeplitz, we take advantage of the Levinson- 
Durbin^ method of solving problems of the form b = Ma, 
where M is a Toeplitz matrix. The Levinson method is 
considered weakly stable, as it is susceptible to numerical 
round-off errors when the matrix is close to degenerate 
(i.e. two or more witness sensors carry nearly identical 
information about the noise source) . For well conditioned 
matricies it is much faster than brute force inversion of 
the matrix^. 

The filtered output of the witness signals is y{n) = 

rp rp 

w^x and e(n) = d{n) — y{n) = d{n) — w^x is the fil- 
tered or minimized target signal, where x is the indepen- 
dent witness signal measuring our noise source. This new 
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minimized signal e represents what the the disturbance 
d would be if we were able to subtract the seismic noise 
before it entered the system. 

Since we are using this feed-forward technique to re- 
duce noise in the length of our cavities, and the cavity 
length carries gravitational wave information, we must 
be careful not to subtract the science signal along with 
the noise. For seismic noise, we are trying to cancel noise 
below 10 Hz, while the gravitational wave signal band 
above 10 Hz. Even so, if we use a filter of the same 
length as the data stream, we could perfectly cancel ev- 
erything, including seismic noise and gravitational wave 
information. However given the short length of our filters 
and the long stretches of data we calculate the subtrac- 
tion over, the distortion of the gravitational wave signals 
should not be significant. 

While the Wiener filtering occurs entirely in the time- 
domain, we examine plots in the frequency-domain of 
the filtered and unfiltered signals to determine the level 
of subtraction achieved. In Section [1111 we describe the 
Wiener simulations on a cavity in our lab. In Section [TVl 
we make similar estimates for one of the 4 km LIGO 
interferometers. Finally, in Section [V| we demonstrate 
the performance of a real-time seismic noise cancellation 
system in the lab. 



III. WIENER FILTERING AT THE 40M 
INTERFEROMETER 

At our 40 m prototype interferometei^ lab at Caltech, 
both static Wiener filtering and adaptive filtering algo- 
rithms have been applied to a suspended Fabry-Perot 
triangular ring cavity's feedback signal. We have used 
2 Giiralp CMG-40T seismometers^ and several Wilcoxon 
731 A accelerometers^^ as our independent witness chan- 
nels (x in Section [ll|, and the low- frequency feedback 
signal for the cavity length as the target channel {d in 
Eq. [T]) to reduce. 

Figure [l] shows the locations of the witness sensors rel- 
ative to the cavity mirrors. The mirrors of the cavity 
are suspended as pendulums with a resonance of ~1 Hz 
to mechanically filter high frequency noise, with the sus- 
pensions sitting on vibration isolation stacks to further 
isolate the optics from ground motion. The 'stacks' are 
a set of 3 legs supporting the optical table on which the 
mirror sits, with each of the legs consisting of alternating 
layers of stainless steel masses and elastomer springs^^EU^ 

We used MatlatP^ to import the data for the length 
feedback signal for our cavity, and to construct and apply 
the Wiener filters. Since the feedback control bandwidth 
is 0:^50 Hz, the feedback signal can be used as an accurate 
measure of the seismic disturbance at low frequencies. In 
Figure |2] we show results of a day-long simulation study. 
This study was done to determine the length of time we 
can use a set of static filters before updating. We use 
1 hour of data to train and calculate a single Wiener 
filter, and then apply that filter to 10 minute segments 




FIG. 1. Locations of seismometers and accelerometers in rela- 
tion to the cavity mirrors. Round trip length of the triangular 
cavity is 27 m. 



of data for one day, using a 31 second long, 2000 tap 
filter with a sample rate of 64 Hz. In Figure |2]3, we 
select a few typical traces to illustrate the capabilities of 
the filter, while in Figure [2^ we show the full results as 
a spectrogram, whitened by normalizing to the spectra 
during the time the filter was being trained. We see large 
amounts of noise reduction both at the broad stack peak 
at ~3 Hz and around the 16 Hz vertical mode of the 
mirror pendula. 

We also include the noise contributions of our seis- 
mometers in Figure [3] to demonstrate how close we are 
able to get to the fundamental limit of Wiener filtering. 
Since the Wiener filter accepts, as inputs, the signals from 
the witness sensors (which have true ground motion in- 
formation plus self-noise of the instruments and noise in 
the readout electronics), all of these noise contributions 
are filtered and added back into our data stream, limiting 
our ability to suppress ground motion below these levels. 
In Figure |3j we show that the differential ground motion 
over the length of the cavity is not much larger than the 
instrument noise of the seismometers. In other words, the 
ground noise over length of the cavity is strongly corre- 
lated below ~ 1 Hz and so the differential motion is much 
smaller than the motion of any individual sensor. Cur- 
rently our measurement of the differential ground motion 
is limited by the apparent instrument noise of the seis- 
mometers, represented by the teal trace in Figure [3| The 
apparent instrument noise is significantly higher than the 
specification, which indicates that there is some unknown 
noise which is uncorrelated between 2 seismometers, even 
when they are placed very close together. We will use 
lower noise sensors and readout electronics and better 
thermal/acoustic isolation of the seismometers in order 
to get better performance on such short baselines. 

The limit to the performance of the feed-forward sub- 
traction seems to be a combination of low frequency noise 
in the seismometers and the feedthrough of noise from 
the auxiliary controls systems of the cavity (e.g. angular 
controls, pendulum damping servos, etc.). 
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FIG. 3. (Color online) Shown are the spectra of the individual 
seismometers (blue dashed and green dash-dot) , the manufac- 
turer's spec for the seismometers' internal noise (purple solid- 
circle), and the differential ground motion along the 13.5 m 
length of the cavity (red solid) . We also show the differential 
noise of the seismometers with the seismometers collocated in 
a stiff seismic vault (teal dash-circle); in principle, this is a 
measurement of the actual seismometer noise floor. It is un- 
known what uncorrelated noise is present in our sensors which 
makes the teal trace so much larger than the specification. 



FIG. 2. (Color online) Result of offline seismic Wiener filter- 
ing on suspended triangular cavity, (a) Spectrogram showing 
the efficacy of a Wiener filter applied offiine over a several 
hour period. Noticeably different traces between ^28 hours 
and ^34 hours are the result of non-stationary anthropogenic 
noise, not a decay of the filter's efficacy, (b) Amplitude spec- 
tral density of the control signal. Dotted red is without sub- 
traction, purple is initial residual, progressively lighter blues 
are 10 hours, 20 hours and 30 hours after filter was trained. 



tests executed on the LIGO interferometers using the 
HEPI actuators during the 6*^ LIGO Science Run will 
be available in a future papei^^. 

Instead of a single cavity, in this case we explored 
the subtraction of seismic noise from the differential arm 
length feedback signal (which is an accurate measure of 
the low frequency ground noise). The sensors are placed 
close to the ends of the interferometer arms and at the 
beamsplitter as shown in Figure [4j 



IV. APPLYING WIENER FILTERING TO A 4 KM LIGO 
INTERFEROMETER 

One of the LIGO sites in Livingston, Lousisiana has 
had a hydraulic external pre-isolator (HEPI) actuation 
system installed since 2004 (the other LIGO site in Han- 
ford, Washington will receive a HEPI system as part of 
the Advanced LIGO upgrade) This HEPI system is 
designed to actuate on the seismic isolation stacks which 
support the suspended LIGO optics to actively reduce 
seismic noise. Initial implementation of the HEPI actu- 
ators only included local seismic isolation between 0.1- 
5 Hz to reduce anthropogenic noise, tidal effects and the 
microseism^^. 

To estimate how the global Wiener filtering technique 
should scale up to a full size interferometer, we analyzed 
data from the 5*^ LIGO Science Ruil^. While this analy- 
sis was done as offline post-processing, results from later 
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FIG. 4. Schematic layout of seismometers relative to interfer- 
ometer mirrors. 
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Figure [5] shows the resulting subtraction efficacy for 
a static filter. The variation in the 0.1-0.3 Hz band 
comes from variation in the ambient level of the dou- 
ble frequency microseismic peak^^. The structure in the 
1-15 Hz band is the usual increase in anthropogenic noise 
during the workday. 

Even though some excess noise is added in the dips 
around 3-5 Hz and 7 Hz, the filter reduces the main con- 
tributors to the RMS of the control signal, and the re- 
duction is remarkably stable over the 30 day timespan. 
This static filter does inject an unacceptable amount of 
noise above 20 Hz, which we will elimiate in the future 
by using more aggresive pre-weighting to disallow such 
noise amplification before calculating the Wiener filter. 

Figure |6] shows the subtraction if we use an acausal 
filter, retraining it every 10 minutes, for the same 30 day 
data set. This filter performs much better than the static 
version. While we cannot apply an acausal filter in re- 
altime, we can utilize causal adaptive filters to achieve 
nearly the same effect as long as the seismic environment 
does not change appreciably on time scales less than 10 
minutes. 

Residuals for both Figures [5] and [6] were calculated us- 
ing 46 second long Wiener filters of 3000 taps at a sample 
rate of 64 Hz. 



V. ONLINE ADAPTIVE FILTERING AT THE 40M 
INTERFEROMETER 

In case the transfer functions between the sensors and 
the target are changing with time, it would be useful to 
use a filter whose coefficients change with time. Such 
an adaptive filter could also take into account changes 
in the 'actuator'. The most simple and common imple- 
mentation of an adaptive filter is the Least Mean Squares 
(LMS) algorithnPSl. 

The Online Adaptive Filtering (OAF) algorithm im- 
plemented at the 40m Lab is the Filtered-x Least Mean 
Squared (FxLMS) algorithnP^. It is based primarily on 
the canonical LMS algorithm; a steepest descent opti- 
mization of a defined error function. Just as in the static 
Wiener filtering in Section [Hj we minimize the RMS of 
the difference between the filtered output and the orig- 
inal feedback signal. The LMS algorithm described in 
Equation [4] takes 'steps' in the direction of the steepest 
gradient until it arrives at a local minima. 



w{n + 1) = w{n) X [1 — r] + /i X e(n) x x{n) (4) 

Here the next iteration's FIR coefficients depend on 
the the current coefficients (w)^ the current witness sig- 
nal (x), the current error signal (difference between the 
target and ffitered signal, e), and the adaptation rate (/i). 
One of the largest challenges with the adaptive filtering 
algorithm is that the success of the algorithm is fairly 
sensitive to the choice of ja. To improve stability against 



■■■■IIIIIIIIHI"^ 

' ■■■■■liwiiRiliriiliilMi'ail 



■ii»''ii'liillifMIPi!il'"WliM 



'ftlr" 


1, 1 1'l 


1 1 III 


iiiyiiiM':'«ih 

1 mil 


iiiiM'f 















I 



I 



5 10 15 20 25 

Time [Days]: tO is GPS 873423660 



M, -10 
^ 10 



% 1 
U 10 

1 

10 

-1 

10 







'*^-^\/-'^ 










Control Signal with no Filtering 


Initial Residual 




10 Days After tO 




20 Days After tO 




30 Days After tO 



10 

Frequency [Hz] 



FIG. 5. (Color online) Result of offline simulated seismic 
Wiener filtering on the 4km LIGO Hanford interferometer, 
(a) Traces are amplitude spectra normalized to the unfiltered 
control signal (red trace in b), which is at a time during the 
filter's training. Filter was trained on 6 hours of data, then 
applied in 10 minute segments. Vertical stripes indicate times 
when the interferometer was not operational. Seismic sub- 
traction is fairly constant on a one month time scale, although 
it is not particularly effective for times when seismic noise is 
significantly different from the training time (b) Selected in- 
dividual spectra from a above. Dotted red trace is before 
subtraction, purple trace is initial residual and progressively 
lighter blues are 10, 20 and 30 days after the filter was trained. 



transients, we modify the usual FxLMS algorithm to in- 
clude a decay constant r. 

The FxLMS algorithm acknowledges that there exist 
phase delays in the path of the target signal which can- 
not be approximated by the LMS method alon^^. To 
account for these phase delays, we filter the incoming 
witness signals with filters identical to those in the tar- 
get signal path. Once we have matched the delays in the 
two different paths, we implement the regular LMS op- 
timization to find the coefficients we will use in our FIR 
filter. The FxLMS algorithm used is sketched in Figure[7| 

We apply the OAF system to the same triangular cav- 



ity as in Section [nil Once again, we use the cavity length 
feedback signal as our targeted signal to minimize, and 
a similar layout of independent witness sensors as shown 
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FIG. 6. (Color online) Result of offline simulated seismic 
Wiener filtering on the 4km LIGO Hanford interferometer, 
using an acausal filter on the same 30 day data set. (a) Traces 
are amplitude spectra normalized to the unfiltered control sig- 
nal (red trace in b). A filter is trained on, and then applied to, 
10 minute segments of data. Seismic noise is more effectively 
suppressed using this constantly updated filter, implying that 
the transfer function is changing on a relatively short time 
scale, and that it is advantageous to update the filter more 
often than once per month, (b) Selected individual spectra 
from a above. Dotted red trace is before subtraction, purple 
trace is initial residual and progressively lighter blues are 10, 
20 and 30 days after beginning. 
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in Figure [T] Unlike Sections III and |IV| which were sim- 
ulations using previously collected data, here we are ac- 
tuating on the cavity in realtime. Figure [8] shows results 
using a 125 second long, 2000 tap filter with /i = 0.01 
and r = 10~^ at a 16 Hz sample rate. The on /oflF traces 
in the adaptive case are similar to estimates made in the 
static Wiener filtering case (Figure [2|. Given enough 
time to adapt, the OAF converges towards the optimal 
filter, but, so far, not completely. Since the adaptive 
system was tested using one Giiralp seismometei^ and 
one Ranger SS-1 seismometei^^, the subtraction is not as 
pronounced as if we had used 2 Giiralps, or other more 
sensitive broadband seismometers. In the next iteration 
of this setup, we will explore the variation in the cancel- 
lation performance as a function of sensor placement. 
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FIG. 7. Block diagram of the FxLMS algorithm used. 



FIG. 8. (Color online) Online Adaptive Filter performance: 
the spectral density of the cavity length fiuctuations are 
shown with the feed-forward on (lower blue trace) and off 
(upper red trace). 



VI. CONCLUSIONS 

We have demonstrated the use of Wiener filter based 
feed-forward seismic noise reduction on a suspended in- 
terferometer. We have also implemented a stable, adap- 
tive feed-forward system which has a performance ap- 
proaching that of the optimal Wiener estimate. These 
techniques can be simply implemented in any general 
laboratory requiring vibration isolation using relatively 
low cost accelerometers and commodity computers and 
DSP software (e.g. Labview). These 'optimal' feed- 
forward schemes operate without having to know a priori 
the transfer function between the disturbance and the 
primary experiment; they can easily be reconfigured to 
adapt to new experimental setups. Similar Wiener filter 
and LMS based techniques have been utilized in other 
experiments, both offline and online, for example isolat- 
ing 2-mirror Fabry-Perot cavities from ground motion"^, 
reducing acoustic noise in oceanography settings^^ and in 
signal processing to decorrelate degenerate witness chan- 
nelPl. 

In the near future, we will work to use this scheme to 
reduce the noise in multiple degrees of freedom of the full 
interferometer. It is clear that this technique can also be 
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applied to remove other sources of environmental noise 
(e.g. acoustic, magnetic, electronic, etc.). For each noise 
source in the gravitational wave band, we will inject soft- 
ware gravitational wave signals into the data stream and 
confirm that they are not distorted. There will certainly 
be new challenges associated with each type of noise, but 
this seems promising as a method which can be employed 
to reduce the influence of environmental noise in a wide 
variety of experimental setups. 
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